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f— *^ ■ Abstract. As well as arising naturally in the study of non-intersecting ran- 

^sj ' dom paths, random spanning trees, and eigenvalues of random matrices, deter- 

minantal point processes (sometimes also called fermionic point processes) are 
relatively easy to simulate and provide a quite broad class of models that ex- 
hibit repulsion between points. The fundamental ingredient used to construct 
a determinantal point process is a kernel giving the pairwise interactions be- 
tween points: the joint distribution of any number of points then has a simple 
expression in terms of determinants of certain matrices defined from this ker- 
nel. In this paper we initiate the study of an analogous class of point processes 

^^ ' that are defined in terms of a kernel giving the interaction between 2M points 

for some integer M. The role of matrices is now played by 2M-dimcnsional 
"hypercubic" arrays, and the determinant is replaced by a suitable general- 

^ <"^ ization of it to such arrays - Cayley's first hypcrdeterminant. We show that 

some of the desirable features of determinantal point processes continue to be 
exhibited by this generalization. 
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j> ■ 1. Introduction 

ly-, ' Motivated by considerations of the behavior of fcrmions in quantum mechanics, 

■^j- , determinantal point processes were introduced in [Mac75 . Surveys of their proper- 

CD ' ties and numerous applications may be found in [DVJ881 ISosOOl |Lyo03| IHKPV061 

^-f.' ! ISTOOl IST03al IST03bl IST04] . 

CD ' We consider a certain extension of this class of point processes. In order to 

motivate our generalization, we first consider a particular case of the determinantal 
point process construction. Suppose that on some measure space (£, .4, /j,) we have 
a kernel K : Y? — > C that defines an L dimensional projection operator for L 2 (/i). 
That is, 
j_j ■ • K(x\ y) = R(y; x) for all x,y £ E, 

• J27j=i K( X H x j) z iZj > for all art, ... , x n & £ and zi,...,z„eC, 

• J" E K(x; y)K(y; z) (i(dy) — K(x; z) for all i,zeE, 

• J s K (x; x) fi(dx) = L. 

The corresponding determinantal point process can then be thought of as an ex- 
changeable random vector with values in S L . The distribution of this random 
vector is a probability measure that has the density 

(si, . . . ,x L ) ~ (LI)- 1 det(K(xf, Xj))t =1 
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with respect to the measure n® L . 

One of the most agreeable things about this construction is that for 1 < N < L 
the JV-dimensional marginal distributions of the random vector have (common) 
density 



(xi, . . . , x N ) ^ (L(L -l)...(L-N+ I))- 1 dek(K( Xi ; x 3 )) 



N 
i,3- 



with respect to the measure /j,® N . Consequently, the conditional distribution of the 
(N + l) st component of the S L -valued random vector given the first N components 
can be computed explicitly (as a constant multiple of a ratio of determinants). It 
is thus possible to simulate the entire E L -valued random vector if one is able to 
simulate a general S-valued random variable from a knowledge of its probability 
density function. 

Various generalizations of determinantal point processes have appeared in the 
literature. Note that 

L 
<r£6[ fe=l 

where &l is the symmetric group of permutations of {1, .... L} and e is the usual 
alternating character on the symmetric group (that is, the sign of a permutation). 
It is natural to replace e by other class functions on the symmetric group (that is, 
by other functions that only depend on the cycle structure of a permutation and 
hence are constant on conjugacy classes of the symmetric group). The most obvious 
choice is to replace e by the trivial character which always takes the value 1, thereby 
turning the determinant into a permanent. Permanental point processes arise in the 
description of bosons and are discussed in [Mac75| fHKP V06| 1DVJ881 IST03a|. IST04] . 
Replacing e by a general irreducible character gives the immanantal point processes 
of [DEOOj . while setting e(a) = a L ~ ly(cr) for -1 < a < 1 and v(a) the number of 
cycles of a gives the alpha-permanental processes introduced in |VJ97j and further 
studied in [ST03J1IHKPV06] . 

All of these constructions have the feature that an exchangeable joint density is 
built up as a linear combination of products of pairwise interactions. In this paper 
we investigate the possibility of building up a tractable joint density as a linear 
combination of products of higher order interactions. In order to accomplish such a 
generalization, it is necessary to have higher order counterparts for both projection 
kernels and determinants. 

Note that K : X 2 — > C is the kernel of an i-dimcnsional projection if and only if 

L 

1=1 

where 4>i , ■ ■ ■ , 4>l are orthonormal in L 2 (/i) . One possible (2M ) th order extension 
of this second order definition is to suppose that: 

• the underlying space S is a Cartesian product Si x • • ■ x Em, 

• the measure /i on £ is a product measure (Xi (8) • • • (g> jj,m, 

• the functions (j> m i : S — * C, 1 < m < M, 1 < £ < L, are given by 
<?W(£i7 ■ ■ • ,%m) = ii>mt{%m), where for 1 < m < M the functions Vw : 
E m — > C, 1 < £ < L, belong to L 2 (fi m ) and are orthonormal in L 2 (fi m ), 
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the kernel K : S 2M — ► C is given by 

L M 

K(y ll . ..,y M ;zi,.. -,z M ) ■= E II ( t ) Tnl{v m )4 > ml{zm) 

1=1 m=l 
L M 



— 2_^ 11 ' t Pml(ymm)lpml(Zmm) ■ 



1=1 m=\ 



Note that the integral 

' M 
|1 4>rnl' (rn){x)(j> m l" (m){x) 



,m=l 
M 



[J,(d; 



X 



Y[ / i J m£'{m){Xm)i>me"(m)(Xm)lJm{dXr, 



is 1 if £'(m) = £"(m) for 1 < m < M, and the integral is otherwise. This 
is analogous to the orthonormality of the functions 4>i, ■ ■ ■ ,4>l appearing in the 
representation above of an L-dimensional projection, and when M = 1 we just 
recover that representation. 

The appropriate generalization of the determinant is given by Cayley's first hy- 
perdeterminant that was introduced in [Cay43| and which we will describe shortly. 
Cayley later introduced other generalizations of the determinant that he also called 
hyperdeterminants and are more natural from the point of view of invariant the- 
ory - see [GKZ92, GKZ94 . Early treatments of the theory related to Cayley's 
original definition may be found in [PasOOl IMui60l IRicl81 IRic301 I01d34cl I01d34bl 
I01d34al I01d361 IOld40| . More recent works are |Sok60l ISok72j . We rem ark that 
Cayley's first hyperdeterminant has been useful in matroid theory |Bar95|, [Gly] and 
we also note the interesting papers [LT03, LT04] in which the calculation of Selberg 
and Aomoto integrals is reduced to the evaluation of hyperdeterminants of suitable 
multi-dimensional arrays. 

Suppose that 

A(ii,. .. ,i M ;ji,-- ■ Jm), 1 < k,- ..,2m, ji,-- -Jm < N. 

is a 2Af-way hypercubic matrix (that is, A is a a 2M -dimensional array or tensor 
that is of the same length, namely N, in each direction). Suppose further that /C is 
a subset of {1, . . . , M}. We define the corresponding hyperdeterminant of A to be 

DetK(A):=^ E E E E IK**)*-*) 



tiSSn &m£&n tiSSjv tmE&n k£K 

N 

n=\ 



JJ_ A((7i(n),. .. ,<TM(n);Ti(ri),.. .r M (n)), 



where &n is the symmetric group of permutations of {1, . . . ,7V} and, as above, 
e is the alternating character. This definition is just the usual definition of the 
hyperdeterminant of a general hypercubic matrix with a general "signancy" , except 
that we have imposed the restriction that the coordinate directions of the matrix 
are grouped in pairs, and each coordinate direction in a pair has the same signancy. 
When M = 1, so that A is just an TV x N matrix, Detjc(A) is either the usual 
determinant or the permanent, depending on whether K is {1} or 0. 
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Note: From now on we will assume that K, is non-empty. 

We are now ready to define a family of exchangeable probability densities. For 
1 < N < L, define the function p^ : Yj N — > C by 

p N ( Xl ,...,x N ):= (QW) M ) Det K (B), 

where B is the 2M-way hypercubic matrix of length N given by 

B(ii, . . . , iulji-, ■ ■ ■ iJm) '■= K(xi 17 . . . , Xi M ; Xj x , . . ., Xj M ). 
That is, 

Pn(xi,. ..,xn) 

'^m^y 1 E ••• E E ■■■ E W^Mr k ) 

N 



II-'*- v^ai (n) :■•■•> ^gm in) i ^n (n) ) ■ ■ ■ ! ^tm (n) J i 



X 

ra=l 

L W f+1 ) E ••• E E ••■ E II **)*>) 



N L M 



II E II ^(^mW^^^TmHm)' 



X 

n=l £=1 m=l 

We will prove the following theorem in Section [3] 

Theorem 1.1. For I < N < L, the function p^ is the density with respect to n® N 
of an exchangeable probability measure on £ . That is, pn > 0, 



Pn{xi,.. .,x n )h® (d(xi,.. .,x N )) = f, 

and pn is a symmetric function of its arguments. The probability measure asso- 
ciated with pn is the common N -dimensional marginal of the probability measure 
associated with pi, . That is, 

p N (xi,.. .,x N ) = / Pl(xi,- ■■ ,x N ,x N+1 ,. .. ,x L ) fi® ( - L ~ N) {d(x N+1 , . . -,x L )). 

The function p^ is the density with respect to n® L of a probability measure on 
£ L = (Ei x • • • x S M ) L ~ £f X • • • X E|^. We show in Section [4] that the marginal 
of this probability measure on £f x • ■ ■ x Ej^, for 1 < M' < M is also given by a 
hyperdeterminantal construction (with the kernel K replaced by a suitable function 
of 2M' variables) as long as K, n {1, . . . , M'} ^ 0. 

If we regard the exchangeable probability measure on E L with density pl as 
the distribution of a point process on E, then it is natural to inquire about the 
distribution of the number of points that fall into a given subset of E. We find 
a relatively simple expression for the factorial moments of such distributions in 
Section 

The key observation behind many of our arguments is an expansion of suitable 
hyperdeterminants that is analogous to the Cauchy-Binet theorem for ordinary 
determinants. This result is an extension of a lemma from Bar95], and we give the 
proof in Section O 



HYPERDETERMINANTAL POINT PROCESSES 5 

2. A HYPERDETERMINANT EXPANSION 

For M = 1, the following result is a consequence of the Cauchy-Binet expansion 
for determinants. (Recall our assumption that /C is non-empty and so our hyperde- 
terminant for M — 1 is a determinant rather than a permanent - the Cauchy-Binet 
expansion for permanents is somewhat different and involves sums over possibly 
repeated indices.) When M > 1 and /C = {1,2, ...,M}, the result is given by 
Lemma 3.3 of }Bar95j . 



Proposition 2.1. Suppose that A is a 2M-way hypercubic matrix with length N 
in each direction that is of the form 



[%l,. -.,lM]Jl,- 



,3m) 



L 



AW(i 1 j)...A( M \i M j)AW(jiJ)---A {M) (JM,£), 



where A^ m ' is an N x L matrix and A^ m ' is the N x L matrix obtained by taking 
the complex conjugates of the entries of A*-" 1 ' . Then Det x. (A) = if L < N and 
otherwise 



Det 



jc\^) = 



E 

c 

E 



J[dct{AP)det(A 



(fch 

c i 



keK 



n i det ( A 



(*0\|2 



.keic 



,(k), 



r(fe) 



J] per (^) Per (^) 

WM2 



I] I Per (40 

k£K 



where the sum is over all subsets C of {1, 2, . . . , L} wzi/i cardinality N and A^ is 
the Nx N sub-matrix of the matrix A^ k > formed by the columns of A^> with indices 
in the set C. 



STEVEN N. EVANS AND ALEX GOTTLIEB 



Proof. We have 
Det,c(A) 

^ j: - j: e ••• e n**M 



AH 



Tk) 



criG&N (TmSSjv ti£Sjv tjjESn fcGAC 

N 

]f[ A(ai(n),.. .,a M (n);Ti(n),. ..r M («)) 



ffi66« ctmG©wti€Sjv t m GSjv fee/C 
V L 

x nE A(1) ( ff i(")^)-" A(M) (^(n),l)i (1) (ri(n),£).-A( M )(rM(n),l) 

n=l £=1 

= m E - E E ••• E n**M7*) 

hi£6n <7m£6n tiGSjv t m SSjv fc£/C 
L L N 

x E"- e n^^W'^)--^^^^)'^) 

£ 1 = 1 l N = \n=\ 

xA^(n(n),e n )---A^(T M (n),£ n ) 

L L 



1= i ^=1 



where 



5(4,4,---,m= E ••• E E ••• E II ^M 7 *) 

(TiGSiv (T M GSjv tiGSiv tm£6« fceAC 



V 



x[]iW(, 1 (n),4)...# f )(, M (T 1 ) 1 4)i (1) (r 1 (r ! ) i g...i( M )(rM(r ! ),g. 



For f = (4, . . . , £jv) G {1, . . . , L} N and m G {1, . . . , M} define an TV x TV matrix 

(m) 



B™ by 



Then 
5(f) = 



BJ.^(i ) j):=^M(i ) ^). 



II ( E e(a)n^ fe )(a(n),^)J ( £ e(r) ]} ^ (fc) (r(»),*n) J 

,fc€AC \o-e6w "=1 / VtG6jv "=1 / 



V 



V 



n e n^v(n)A))( e n^ (fc) (rN,4 

fc^AC \<re©jv "=1 



VtGSjv n=l 



ndet(Bf)det(Bf) 



.fee ac 



n P er(BfVr(Bf) 

J Lfc^AC 



nidet(B}' 



(fc)\|2 



.feeAC 



Hi per (B<' 



(fc)M2 
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Note that the rightmost product is zero unless the entries of the vector £ are 
distinct, because in that case each of the matrices B- for k £ /C will have two equal 
columns and hence have zero determinant (recall that K, is non-empty). Moreover, 
if £'=(£[,.. .,£' N ) and £" = (£'{, ... ,£'^) are two vectors with distinct entries such 
that {£[,..., £' N } = {£'{, ...,l%} = £, then 



|det(B} fcJ )| 2 = |det(^ 

and 

|per(4 fc) )| 2 = |per(4 fc) )| 2 

for all k, because permuting the columns of a matrix leaves the permanent un- 
changed and either leaves the determinant unchanged or alters its sign. 

The result now follows, because for any subset £ of {1, 2, ... , L} with cardinality 
N there are N\ vectors £= (4, . . . , £ N ) with {£i, . . . ,£ N } = £. □ 

3. Proof of Theorem 11.11 
By definition, 

] ^J(Nl) M PN (x 1 ,...,x N )= Det K (B), 

where B is the 2M-way hypercubic matrix of length N given by 
l(ii, . . . ,i M \ ji, ■ ■ -Jm) = K(x ll} . . .,x lM ;x h ,. . -,x JM ) 

L M 

1=1 m=l 
L 

= Y / B {1 HH,£)---B^(t M ,£)B^(j 1 ,£)---B^\j M J), 
e.=\ 

and the N x L matrix i?( m ) is given by 

By Proposition 12. 1[ 
' L 



v AN\) M p N {xi,...,x N ) = J2 

C Lfee/C J Lk4K. 



nidet(Bf)| 



(fe)\|2 



niper(4 fc) )| 



(fc)\|2 



where the sum is over all subsets £ of {1, 2, ... , L} with cardinality N and B,. 
is the N x N sub-matrix of the matrix B^> formed by the columns of B^ with 
indices in the set C. 

It follows that pn(xi, ■ ■ ■ ,Xn) > 0. Also, since the value of the permanent a 
matrix is unchanged by a permutation of the rows and the value of a determinant 
is either unchanged or merely changes sign, the function p^ is unchanged by a 
permutation of its arguments. 
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We have 



{N\) M p N {x 1 ,...,x N ) 



E 



A' 



\\ EE e ( (7 ' i; ) e ( Tfe ) II ( t ) ka k {n){Xn)<t>kT k (n){x n ) 
.keK. <Jk i~ k n=l 



N 



II EE II 0k<T*(n)(Zn)0fcT*(n)(z„) 

Lk^K. cr k r fc n=l 

M N 

= EE"'EE'"EII e KM T fc) II '\\.^m<7 m (n){x n )4> m r m (n){x n ), 

C cri um T\ tm keIC m=ln=l 

where o^ and Tfc in the summations range over bijective maps from {1, . . . , TV} to 
C and e is interpreted in the usual way for such a bijection. 
Now the integral 

' M 

J J_ < Pm(T m (n)\Xn)'PmT m (n){%n) 
,?n— 1 
M 



fl(dx n 



m— 1 



J^J / Vmo m (n) \Xnm)VmT m (n)\^nm) l^\"'X nrn ) 



is equal to 1 if and only if o~ m (n) = r m (n) for 1 < m < M, and otherwise the 
integral is 0. Hence 



E- v 



PAr(a;i,. .. ,x N )/j,® N (d(xi,.. . ,xjv)) = 1 



and 



as required. 



p N (xi, . . . , x N ) n{dx N ) = Pn-i(xi, ..., x N -i), 



4. Varying the order M 

Beginning with a suitable kernel K : S 2M — (Si x • • • x S M ) 2M — > C, we have 
a built a family of functions p^ 1 1 < N < L, where pn is a probability density on 
(Si x • • • x Sm)" ~ Sf x • • • x Ej^ with respect to the measure (<g>™ =1 /i m ) 0N ~ 
m=1 Hm N - F° r 1 ^ -^' < ^j ^ is natural to ask about the push- forward of the 
probability measure corresponding to the density pn by the projection map from 
0^1/4^0®^^ given by 

{Xnm)l<n<N, l<m<M *-* {x n m)l<n<N, \<m<M' ■ 

The answer is given by repeated applications of the following result. 

Theorem 4.1. Suppose that either M $ K, or M e K, and K\ {M} ^ 0. Set t := 

K \ {M}, t ■- n^ 1 S m , and jX = (g)^ 1 Mm- £>e/?ne a kernel K : ± 2 ( M -^ -> C 

L Af-l 

K(Vl, ■ ■ ■ ,VM-i;Zl, ■ ■ ■ ,ZM-l) : = E II < i ) ™-(-(yrn)(i ) mt{Zm)- 

1=1 m=l 
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The function 

(x n m)l<n<N, l<m<M-l h ^ PN{{%nm)l<n<N, l<m<M-l) 

PN{{ x nm)l<n<N,l<m<M) Hm (^( X 1M , ■ ■ ■ , Xnm)) 



is a probability density with respect to the measure fj,® N . The probability density pN 
is constructed from the kernel K and the set of indices K, in the same manner that 
the probability density p^ is constructed from the kernel K and the set of indices 
K. 

Proof. As in the proof of Theorem ll.l[ 



.{NX) PN{{Xnm)l<n<N,l<m<M) 

M N 

= 5Z5Z'"^zx]'"X! n e (°'*) e ( 7 *) n n ^^(n)^™)^^^)^™), 

C a i a M Tx tm keK. m=ln=l 

where o~k and Tk in the summations range over bijective maps from {1, . . . , TV} to 
C. 

Observe that the integral 



jv 



] _[ i>Mcr M (n) {XnM)i>MT M (n) {%nM ) 



Hjf {d{xiM,-- -,xnm)) 



is 1 if and only if Cm (n) = Tif(n) for 1 < n < N (that is, if and only if Cm — Tm), 
and the integral is otherwise. For each choice of the set £, there are N\ choices 
of the pair of bijections (<7mt t m) such that o~m = Tm, and for all of these choices 
we have, of course, that c(<7m) = c{tm) if M G £. 
It follows that 

( n)( N '^ M / PN{{Xnrn)l<7i<N,l<7n<M) [J.M N (d(x 1M , ■ ■ .,X N Al)) 

C d\ (TM — 1 T\ TM-1 k£K 

M-l N 

x - J[ ^[ j[ J^ V- , ma rn (n) \Xnm)V^mr m (n) [Xnm )■, 
m—1 n— 1 

as claimed. □ 



5. The number of points falling in a set 

Write (_Xi, . . . ,Xl) for a £ L -valued random variable that has the distribution 
possessing density pl with respect to the measure /i® L . Given a set C £ A, let 
Jc ■= #{1 < £ < L : Xg E C}, so that Jc is a random variable taking values in 
the set {0, 1, . . . , L}. The distribution of the random variable Jc is determined by 
the factorial moments E[J C (J C - 1) ■ ■ ■ {J C - N + 1)], 1 < N < L. 
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If we write It for the indicator random variable of the event {X( G C}, then 
Jc = h + ■ ■ ■ + II- By the exchangeability of (Xi, . . . , Xl), we have 

E[J C (J C - 1) ■ ■ ■ (J c - N + 1)] 



,x N )). 



= L(L- 


_!)... 


(L 


-N + l) 


E[Ji 


■■In] 






= L(L- 


_!)... 


(L 


-N + l) 


/ . 

Jc N 


D N (X 1 ,... 


,X N )fl C 


m (d( Xl , 


Suppose now 


that C 




C x x • • • 


x C M 


, with C m 


c y 

_ 't-'rm 


1 < m < 


E[J c (Jc - 1) 


■■■(Jc 




N + l)] 











M N 



= (N\)- M+ 



x ee ee En««nn^, 



(m) 



i(n),r rn (n) ' 



where i/( m ) is the L x L matrix defined by 



H^ m \£J') 



ip m £(y)4>m£> (y) Hm(dy) , !<£,£' <L. 



Thus, 

nj c (Jc 



1)---(J C -N + 1)]=N\J2 



n det (^ fc) ) 



£ LfceK 



[] per (iff) 



where H c is the N x N sub-matrix of H ( k > with rows and columns indexed by C. 

Remark 5.1. When M = 1 (so we are dealing with a determinantal point process), 
the last expression is just the trace of the ( N ) x ( N ) compound matrix consisting 
of the minors of H^ 1 ' with N rows and columns. This compound matrix has as its 
eigenvalues all the products of the eigenvalues of H^ 1 ' taken N at a time, and so its 
trace is the sum of all such products. This shows that Jc is distributed as the sum 
of L independent Bernoulli random variables that have the eigenvalues of H^> as 
their respective success probabilities - a result that appears in |ST03bl [HKP V06J . 
We have been unable to find an analogous probabilistic representation for Jc for 
general M. 
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